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Abstract 



The 7-ray light curves and spectra are presented for a set of theoretical Type la 
supernova models including deflagration, detonation, delayed detonation, and 
pulsating delayed detonations of Chandrasekhar mass white dwarfs as well as 
merger scenarios that may involve more than the Chandrasekhar mass and he- 
lium detonations of sub- Chandrasekhar mass white dwarfs. The results have 
been obtained with a Monte Carlo radiation transport scheme which takes into 
account all relevant 7-transitions and interaction processes. The result is a set 
of accurate line proflles which are characteristic of the initial ^^Ni-mass distri- 
bution of the supernova models. The 7-rays probe the isotopic rather than just 
the elemental distribution of the radioactive elements in the ejecta. Details of 
the line proflles including the line width, shift with respect to the rest frame, and 
line ratios are discussed. With sufficient energy and temporal resolution, dif- 
ferent model scenarios can clearly be distinguished. Observational strategies are 
discussed for current and immediately upcoming generations of satellites (CGRO 
and INTEGRAL) as well as projected future missions including concepts such 
as Laue telescopes. With CGRO, it is currently possible with sufficiently early 
observations (near optical maximum) to distinguish helium detonations from ex- 
plosions of Chandrasekhar mass progenitors and of those involving mergers up to 
a distance of about 15 Mpc. This translates into one target of opportunity every 
eight years. SNe la up to about 10 Mpc would allow detailed CGRO studies of 

line ratios of ^^Co lines. INTEGRAL will be able to perform detailed studies 
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of the Co line profiles with a range comparable to CGRO. The superior sensi- 
tivity of INTEGRAL for low energies makes detection and detailed study of the 
positron annhilation line and appropriate low energy ^^Ni lines possible up to 
about 10 - 15 Mpc for all models. This capability means that this lower energy 
range may be the most useful for INTEGRAL detection and study of SNe la. 
Such studies will allow the determination of the precise time of the explosion. 
Whereas the current generation of 7-ray detectors will allow the study of super- 
novae which are discovered by other means, a new generation of proposed 7-ray 
detectors with sensitivity of about 10 photons sec cm would generate 



the opportunity to discover supernovae by their 7-ray emission up to a distance 
of ~ 100 Mpc. This would aUow a systematic study of the variety of SNe la in 
terms of their 7-ray properties, independent of their optical properties. In addi- 
tion, since 7-rays are not obscured by the host galaxy, such experiments would, 
for the first time, provide absolute supernova rates. Relative rates as a function 
of the morphology of and position in the host galaxy could be studied directly. 

Subject headings: Supernovae and supernovae remnants: general - gamma rays 

1. Introduction 

Type la Supernovae (SNe la) are among the most spectacular explosive events. 
They are major contributors to the production of heavy elements and hence a 
critical component for understanding the life cycles of matter in the Universe 
and the chemical evolution of galaxies. While much progress has been made on 
understanding the underlying physics of the thermonuclear explosion (Khokhlov, 
Oran and Wheeler 1997ab; Niemeyer & Woosley 1997), the basic processes re- 
quire further elucidation and testing and the progenitor evolution is still a deep 
mystery (Wheeler 1996). The great brightness of SNe la has made them a valu- 
able tool for the measurement of extragalactic distances (Freedman et al. 1994) 
and to determine the shape of the Universe (Perlmutter et al. 1995, Schmidt et 
al. 1997). The absolute brightness must be known either by using other distance 
calibrators (Hamuy et al. 1996, Ries, Press, & Kirshner 1996, Sandage & Tam- 
mann 1996) or theoretical models for the light curves and spectra (Hoflich 1995, 
Hoflich & Khokhlov 1996, HK96 hereafter). A detailed physical understanding of 
SNe la is important for its own sake. In addition, such a physical understanding 
is necessary to constrain systematic errors when using them at high rcdshifts 
to determine the deceleration parameter qq and other cosmological parameters. 
Evolutionary effects on the progenitor population, for instance the initial IMF 
and the metallicity, are expected to produce systematic errors (Hoflich et al. 
1997). 

It is widely accepted that SNe la are thermonuclear explosions of carbon- 
oxygen white dwarfs; however, rather different progenitor models are under dis- 
cussion. Three main scenarios can be distinguished. The first group consists of 
carbon-oxygen white dwarfs with a mass close to the Chandrasekhar mass which 
accrete mass through Roche-lobe overfiow from an an unevolved (CV-like) or 
evolved companion star. In these accretion models, the explosion is triggered 
by compressional heating. From the theoretical standpoint, the key questions 
arc the condition under which the thermonuclear explosion starts and how the 
flame propagates through the white dwarf. The ignition and propagation may be 
affected by the rotation and mass accretion rate of the white dwarf which will in 
turn be affected by the metallicity, age, and population of the progenitor system. 



The second group of progenitor models consists of two low-mass white dwarfs in 
a close orbit which decays due to the emission of gravitational radiation. This 
eventually leads to the merging of the two white dwarfs. In this scenario, the to- 
tal mass of the merged object may exceed the Chandrasekhar mass and rotation 
may play an especially significant role. A third class of models are those involv- 
ing helium detonations, i.e. double detonation of a C/0- white dwarf triggered 
by the detonation of an outer helium layer in low- mass white dwarfs. 

Each of these scenarios will have different implications for the resulting super- 
nova statistics, the evolutionary ages of the progenitor systems and for chemical 
evolution of galaxies and the Universe. 

Gamma rays are of particular interest as a diagnostic tool of these various 
progenitor models because they allow the direct observation of radioactive iso- 
topes which power the observable light curves and spectra. Almost three decades 
ago Clayton et al. (1969) pointed out that the 7-ray lines from SN la explosions 
should be detectable at distances up to a few Mpc. Since then many authors have 
studied the gamma-ray signature of SN la (Clayton 1974; Ambwani & Sutherland 
1988; Chan & Lingenfelter 1988, 1990, 1991; Burrows & The 1990; Burrows et 
al. 1991; Hoflich, Khokhlov & Miiller 1992, 1993ab, 1994, Kumagai & Nomoto 
1997). These studies have shown that 7-ray observations of SNe la provide a 
powerful diagnostic tool to test proposed theoretical models. The key role of 
7-rays for understanding of SNe la became obvious in the case of the peculiar 
SN 1991T. Because of the lack of 7-ray detection and/or firm upper limits, the 
nature of SN1991T was not fully revealed despite various efforts (Hoflich et al., 
1994, Kumagai & Nomoto 1997) and an excellent set of data from the optical to 
the IR (e.g. Meikle et al. 1996) and even a distance from 5-Cephei stars (Saha et 
al. 1996). Nevertheless, this event clearly showed the potential and implications 
for 7-ray observations (e.g. Lichti et al. 1994; Leising et al. 1995, and references 
therein) . 

Although systematic studies of the 7-ray emission properties of different sce- 
narios have been published, several scenarios, some of them with rather strong 
observational and physical support, have not been considered (e.g. Burrows & 
The 1990; Chan & Lingenfelter 1990, 1991; Hoflich et al. 1992) or where both 
Chandrasekar and sub-Chandrasekar type models have been included, no detailed 
study of line profiles has been performed (e.g. Kumagai & Nomoto 1997). 

In addition, little attention has been paid to specific observational strategies 
and the information required to allow instrument designers to estimate the per- 
formance of their instruments and to analyse future data. A detailed knowledge 
of the line position is critical for instruments which provide only a very limited 
wavelength range such as the Laue telescopes currently in the design phase. Even 
for conventional 7-ray telescopes such as CGRO or the upcoming INTEGRAL, 
the detection limits may improve if information on the line position and profile 
is taken into account. Due to the low signal to noise ratio, the detectibility de- 



pends sensitively on the line profile and on the knowledge of the line position. 
In principle, the information from the models and the measurement can be com- 
bined during the data analysis to filter out the high noise rather than comparing 
predicted and observed fiuxes after the reductions. 

This paper will seek to fill some of the gaps just mentioned. In the next section 
we describe briefly the numerical methods and the hydrodynamical models on 
which this study is based. In § 3 the properties of line profiles are quantitatively 
discussed along with the 7-ray light curves. In § 4, the expected rate and the 
possible applications for supernova searches with current and future instruments 
are briefly discussed. § 5 presents a flnal discussion and conclusions. 



2. Brief Description of the Numerical Methods 

2.1. Hydrodynamics 

The dynamical explosion models are calculated using a one-dimensional La- 
grangian hydro code with an artificial viscosity (Khokhlov, 1991) and a radiation- 
hydro code that includes nuclear networks (HK96 and references therein). The 
latter code is based on the light curve code presented in (Hofiich, Miiller & 
Khokhlov 1993a) that solves the hydrodynamical equations explicitly by the 
piecewise parabolic method (CoUela and Woodward 1984). This code includes 
the solution of the radiation transport implicitly via the moment equations, ex- 
pansion opacities, and a detailed equation of state. Typically, 300 to 500 depth 
points are used. The 7-ray transport is omitted during the hydrodynamical phase 
because of the high optical depth of the inner layers. Nuclear burning is taken 
into account using Thielemann's network (Thielemann, Nomoto & Hashimoto 
1994 and references therein). During the hydro, an a-network of 20 isotopes is 
considered to properly describe the energy release. The final chemical structure 
is calculated by postprocessing the hydrodynamical model using a network of 
216 isotopes. 

2.2. Gamma-ray transport 

For a given explosion model, 7-rays are calculated using a detailed Monte Carlo 
7-ray transport scheme (Hoflich et al. 1992, 1993b). A Monte Carlo scheme has 
been used because of its flexibility with respect to the treatment of physical 
processes, because relativistic effects can be included in a straightforward way, 
and because it allows for a consistent treatment of both the continua and lines. 

Photons arc initialized in the transport code according to the distribution of 
^^Ni and ^^Co in the hydrodynamical model. All ^^Ni (6) and ^^Co (46) de- 
cay lines have been taken into account. For the (3^ decay, the assumption is 
made that all positrons form positronium at the point where they are produced 



by radioactive decay and hence that they decay locally to produce annihilation 
photons. The positron leakage can be neglected during the first few months, but 
it may become important during late stages, especially for low mass envelopes 
or if the Ni is located in the outer layers. For a discussion of positron leakage, 
see Colgate, Fryer & Hand (1997) and references therein. About 1 — 5 10^ initial 
decay photons are created for each Monte Carlo calculation. The direction of 
the photons is chosen randomly in the comoving frame of the expanding mat- 
ter. Following Ambwani & Sutherland (1988), we have included the following 
three interaction processes in the Monte Carlo transport code: Compton scatter- 
ing according to the full angle-dependent Klein-Nishina formula, pair-production 
with cross section taken from Hubbell (1969), and photoelectric absorption with 
cross sections determined on the basis of the tables given by Veigle (1973). The 
emergent 7-ray spectra have been calculated with a high energy resolution of 
X/SX — 600. This resolution is sufficient to obtain accurate line profiles and 
fluxes. For more details see Hoflich et al. (1992). 



3. Dynamical Models 

The 7-ray spectra have been computed for a large variety of models including 
delayed detonations (the DD-series of Hoflich et al. 1997), pulsating delayed 
detonations (the PDD-series presented in Khokhlov 1992, Hoflich, Khokhlov & 
Wheeler, 1995, and HK96), binary white dwarf merger scenarios (the DET2env 
series of Khokhlov, Miiller & Hoflich 1993), and helium-detonation models (the 
HcD scries of HK96). The classical detonation scenarios are omitted because this 
class of models is known to produce hardly any intermediate mass elements and 
hence to fail by inspection the criterion of agreement with spectral observations. 

The flrst group of models we consider (the delayed detonation and pulsat- 
ing delayed detonation models) consists of carbon-oxygen white dwarfs with a 
mass close to the Chandrasekhar limit. These systems presumably result from 
accretion through Roche-lobe overflow from a companion star (Nomoto & Sugi- 
moto 1977; Nomoto 1982). In these accretion models, the explosion is triggered 
by compressional heating and the heat wave traveling inwards from any nuclear 
shell-burning region at the surface. Typically, the thermonuclear runaway occurs 
at central densities pc ~ 2 — 3 x 10 g cm . In the innermost regions, burning 
occurs at low Ye and, consequently, neutron-rich isotopes are produced rather 
than ^^Ni as shown in Fig. 1 (panel 1). 

The second group of progenitor models is based on the scenario of two low-mass 
white dwarfs in a close orbit which decays due to the emission of gravitational 
radiation thus leading eventually to the merging of the two white dwarfs as the 
smaller mass, larger radius white dwarf fills its Roche lobe first and undergoes un- 
stable mass transfer (Iben and Tutukov 1984, Webbink 1984, Iben and Tutukov 
1997). The combined mass of the two white dwarfs exceeds the Chandrasekhar 




Figure 1: Some abundances and Ni distribution as a function of mass in stellar 
units for the delayed detonation model DD200c, the merger model DET2env4 
and the helium detonation model HeDlO. 



limit, but the larger mass white dwarf which receives the mass is nevertheless 
itself less than the Chandrasekhar limit. Its subsequent carbon ignition is ex- 
pected to occur at low density as mass accretes from the extended C/0 envelope 
formed by the disrupted star (Benz et al. 1989). Since burning occurs at low 
densities, there is little neutronization and ^^Ni is produced in the center (Fig. 
1, panel 2). 

The third class of models explored here are based on low mass white dwarfs 
in which a degenerate helium layer builds up on top of a degenerate C/0 core. 
The explosion is triggered by detonation of the helium layer which then, in some 
circumstances, causes the inner C/0 core to detonate as well. This scenario 
for normally bright SNe la was explored by Nomoto (1980), Woosley, Weaver & 
Taam (1980), and most recently by Livne & Glasner (1990), Woosley & Weaver 
(1994), Livne & Arnett (1995), and HK96. This scenario has also been suggested 
as an explanation of the subluminous SNe la (Woosley & Weaver 1994, Ruiz- 
Lapuente et al. 1993). In this class of models, the thermonuclear runaway 
starts near the bottom of the outer He layer. The Ni distribution, and the 
subsequent 7-ray emission differs significantly from the models based on the 
other scenarios. Due to the large energy release during He burning, ^ O.IMq of 

Ni is produced in the outer region originally composed of helium. Subsequently, 
a detonation travels inwards and ignites the C/0 core at low density, yielding a 



central concentration of ^"Ni at lower velocity. This bi-modal "^"Ni distributions 
differs significantly from the that expected in the other scenarios (Fig.l, panel 3). 
Depending on model parameters, the amount of this inner Ni can be substantial 
(see model HeDlO in Table 1) and hence represent a potential normally bright 
SN la or it can be as small as the outer layer, O.IMq (see model HeD6 in 
Table 1), which could represent a subluminous event. 



TABLE 1 



Parameters of Dynamical Models of SNe la. 



Model 


Mode of 




Pc 


a 


Ptr 


Mm 




explosion 


[Mq] 


[10^] 




[10^] 


[mq] 




DD200C 


delayed det. 


1.4 


2.0 


0.03 


2.0 


0.613 


DD201C 


delayed dot. 


1.4 


2.0 


0.03 


1.7 


0.56 


DD202C 


delayed det. 


1.4 


2.0 


0.03 


2.5 


0.67 


DD203C (2:3) delayed det. 


1.4 


2.0 


0.03 


2.0 


0.59 


PDD5 


pul.del.det. 


1.4 


2.7 


0.03 


0.76 


0.12 


PDD6 


pul.del.det. 


1.4 


2.7 


0.03 


2.2 


0.56 


HcD6 


Hc-dct. 


0.6+0.172 


.0091 






0.252 (.08) 


HcDlO 


He-dct. 


0.8+0.22 


.036 






0.75 (0.1) 


DET2env4 


det.+envelope 


1.2 + 0.4 


0.04 






0.63 



Basic quantities of the dynamical models are given in Table 1. The quantities 
given in columns 3 to 7 are: M-k, white dwarf mass; pc, central density (in CGS) 
of the white dwarf; ct, ratio of deflagration velocity and local sound speed; p^^, 
transition density (in CGS) at which the deflagration is assumed to turn into 
a detonation; M^^, mass of synthesized ^^Ni . The number in brackets is the 
amount of high velocity Ni (Fig. 1). For the helium detonations and the 
envelope (merger) models, the mass of the C/0 core, and the He-layers of the 
hydrostatic white dwarf or of the extended CO-envclope, respectively, are given 
separately. For all models the initial mass ratio of C/0 was assumed to be 1 
except for model DD203c, for which C/0 was assumed to be 2/3. 

The models for delayed detonations, pulsating delayed detonations and enve- 
lope (merger) models are selected from the total array of those models to be 
the subset with parameters which produce reasonable multi-color light curves of 
events in the observed sample of both normally bright and subluminous SNe la. 
Calculations of spectral evolution or even maximum light spectra do not yet exist 
for the whole sample of models, but the selected models are chosen with com- 
position/velocity profiles that are not in obvious disagreement with observations 



(as arc the pure detonation models). For comparison, we note that the structure 
of the deflagration model W7 (Nomoto et al. 1984) is similar to that of DD201c. 
The helium-detonation models have been advocated in the literature (Livne & 
Arnett 1990, Woosley and Weaver 1994, Arnett 1997) as an explanation for both 
normally bright and subluminous SNe la. Sub-Chandrasekhar progenitors would, 
in principle, be easier to produce from stellar evolution than Chandrasekhar-mass 
progenitors. There are problems with both the light curves and spectra of this 
class of models. The light curves are generally too blue (Hoflich et al. 1997), and 
spectra are dominated by Co/Ni. In addition, the normally bright models show 
excessively weak Si II at maximum light and the subluminous versions of the 
model have the Si moving at too low an expansion velocity (Hoflich et al. 1997, 
Nugent et al. 1997). Despite these objections, it is important to investigate this 
class of models and subject it to observational check in as many ways as possible. 
As we will argue below, future 7-ray observations can provide a critical check of 
this class of models. 



4. Spectral evolution 

The evolution of the hard X-ray and 7-ray spectrum of the delayed detonation 
model DD202c is shown in Figs. 2 and 3 (note the changing scale on the vertical 
axis). The spectrum and its evolution is dominated by nuclear lines and the 
decay ^^Ni — (j^ _^56 respectively. During the flrst 2 weeks, ^^Ni lines 
dominate. The continuum is formed by Compton scattering produced by "down" 
scattered line photons. Below 0.3 MeV (Fig. 3), the spectrum is dominated by 
the Compton scattering continuum with the exception of a very strong Ni line 
at 0.16 MeV which can be seen for at least 40 days after the explosion. This line 
may be observable with a new generation of hard X-ray instruments (see section 
5). 

After two weeks or so, the ^^Co lines become the most prominent line features. 
With increasing time, the Compton optical depth decreases and, eventually, the 
Compton continuum becomes weak (lower panels in Fig. 2). Due to the two- and 
three-photon decay of ortho- and para- positronium, a strong 511 keV line and 
a flat continuum below 511 keV is produced. After about 6 weeks, the overall 
spectral distribution remains unchanged. 

The spectral evolution of other evolutionary and dynamical scenarios is shown 
in Fig. 4 which illustrates model PDD6, a pulsating delayed detonation model 
which is an alternative for a normally bright SNe la and HeDlO, which is an 
example of a bright helium shell, double-detonation model. Qualitatively, the 
profiles look rather similar for all models (including deflagrations like W7, Hoflich 
et al. 1992), with the noticeable exception of the helium detonations. For the 
Chandrasekhar mass or merger models, the model spectra differ mainly at late 
times when the ejecta are transparent. At this phase, the central hole in the 
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Figure 2: Evolution of the gamma-ray spectra of the delayed detonation model 
DD202c. The dashed lines corresponds to the contribution of ^^Ni. The strongest 
lines are ^^Ni (0.48MeV), the positron decay line at O.SlMeV, ^^Ni (O.SlMeV), 
^^Co (0.84MeV) and ^^Co (1.23MeV). 




Figure 3: Same as Fig. 2, but for lower energies. The strongest lines are the 
^^Ni lines at 0.158, 0.26 and 0.48 MeV and the positron decay line at 0.51 MeV. 

^^Ni distribution results in flat line profiles. Note the difference between model 
DD202c in Fig. 2 and model PDD6 in Fig. 4 at 231 days, where the latter shows 
a distinct flat top to the lines. These differences give promise for constraints on 
the central density of the progenitor, the initial deflagration speed and, possibly, 
mixing processes (Hoflich et al. 1992). In the absence of gravitation at the center, 
the flame speed is expected to be close to the laminar deflagration speed and then 
to accelerate as Rayleigh- Taylor instabilities take over (Khokhlov 1995, Khokhlov 
et al. 1997ab, Niemeyer & Hillebrandt 1995, Niemeyer & Woosley 1997). Late 
time 7-ray spectra provide valuable information on the velocity distribution of 
radioactive material in the expanding envelope and, thus, they probe the central 
density of the WD. To differentiate between different Chandrasekhar mass and 
merger models requires high resolution spectra corresponding to about 1000 km 
sec"-'^ because the maximum ^^Ni velocity varies from about 2000 to 3000 km 
for models which are able to reproduce observations. The ^^Ni abundance 
typically changes by about a factor of 80 percent over a velocity range of 1000 
km s~^. Thus, signal to noise ratios of the order of 5 are required to separate 
the various effects. An alternative, less observationally demanding method is 
discussed below. 

Helium detonations show signiflcantly different 7-ray spectra and spectral 
evolution compared to the Chandrasekhar mass and merger models. At all 
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Figure 4: Same as Fig. 2, but for the normal bright supernova models PDD6 
and HeDlO 

phases, the helium-detonation models can thus be distinguished from both Chan- 
drasekhar mass and merger models with appropriately sensitive observations. 
The Ni at high velocities from the outer detonated helium shell results in very 



broad profiles beginning a few days after the explosion. Two to three weeks after 
explosion, near or after maximum optical light, the escape probability of 7-rays 
from the central, low velocity ^^Ni region is sufficient to produce a narrow com- 
ponent atop the broad component (see panel 2 for model HeDlO at 23 days in Fig. 
4). Depending on the total amount of ^^Ni , the central component dominates 
after about one month and, eventually, a narrow strongly peaked line remains 
with low amplitude, but very extended wings. At late times, the amplitude of 
the high velocity component is lower than the low velocity component by almost 
an order of magnitude even though the ^^Ni mass is lower by a factor of only 
about 7 in the normally bright model (HeDlO) and 3 in the subluminous model 
(HeD6). The reason is that the the outer ^^Ni is smeared out over about twice 
the velocity range, thus diluting the flux at a given energy. 

The differences between the Chandrasekhar mass and sub-Chandrasekhar mass 
models are even more defined when the subluminous versions of each are com- 
pared. In massive white dwarf models, the mean expansion velocity of the 
^^Ni decreases with the decrease in total ^^Ni mass. This means that the sublu- 
minous models that have less nickel according to the dynamical models (Hoflich, 
Khokhlov, and Wheeler 1995) have narrower lines compared to the models for the 
normally bright SNe la (see also Fig. 7). For sub-luminous helium detonations 
at early epochs, the ^^Ni at high velocities makes the dominant contribution to 
the flux in a given ^*^Ni line. Consequently, the ^^Ni lines are very broad and 
nearby lines are blended, unlike the Chandrasekhar-mass models. The mean ve- 
locity of the inner ^^Ni in the brighter helium-detonation model (HeDlO in Fig. 
4) is higher than that for the less luminous model (HeD6 in Fig. 5). The velocity 
profiles of the outer helium layer are similar for the two models. As a result, at 
late times, the broad "shoulder" remains more visible in the 7-ray spectra of the 
subluminous helium detonation models, as can be seen by comparing the spectra 
at 231 days of model HeDlO in Fig. 4 with model HeD6 in Fig. 5. This "shoul- 
der" in the subluminous models is in especially strong contrast to the late-time 
line profile of the subluminous Chandrasekhar-mass model as shown in Fig. 5. 

In principle, the outer He- Ni layers of helium-detonation models may also 
be observable in the optical and IR for the first few months, because non-thermal 
electrons can provide excitation of the corresponding lines (Hofiich et al. 1997). 
In the absence of a magnetic field, the positrons can freely escape after a few 
months and the optical and IR emission would decline. With a tangled mag- 
netic field and a continued trapping of the positrons, this associatedl optical 
and IR emission could continue more strongly to later epochs. Early optical and 
IR spectra may thus provide an additional means to detect the He- Ni layer, 
but this possibility is currently controversial. Whereas two groups (Hofiich et 
al. 1997, Nugent et al. 1997) find that helium detonations give maximum light 
spectra and light curves that do not reproduce the observations of SNe la. Pinto 
& Eastman (1995, private communication and conference talk in Spain) argue 
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Figure 5: Same as Fig. 2, but for the subluminous supernovae models PDD5 
and HeD6. 

that agreement is possible. Consequently, 7-rays with their straightforward in- 
terpretation must be regarded as the best way to ultimately clarify the nature 
of SNe la. 



^ 10 



- 



' I ' ' ' ' I 
DD201C 




_l I I I I L 



100 200 
time [days] 



10 



= 



' I ' ' ' ' I 
DD200C 




100 200 
time [days] 




100 200 
time [days] 



^ 10 - 




100 200 
time [days] 



^ 10 



5 - 



0) 

c 

- 



' I ' ' ' ' I 
HeDlO 




100 200 
time [days] 




100 200 
time [days] 




100 200 
time [days] 




100 200 
time [days] 



_0.84 MeV ^ Co 



.1.24 MEV ^ CO 



Figure 6: Time evolution of tlie line shift relatively to the rest frame. Note that 
the line shift is insensitive to the specific line considered. 



4.1- Quantitative Analysis of Line Profiles 

In order to provide specific guidance for observational missions, it is appro- 
priate to present more specific detail of the properties of the 7-ray line profiles 
given in the previous section. Of special importance are the line widths and 
the the displacement of the line with respect to the rest wavelength due to the 
optical depth and structural effects in the models. The time evolution of the 
displacement of the line for a variety of models is given in Fig. 6 for the two 
strong lines of ^^Co at 0.84 and 1.24 MeV. Fig. 7 shows the evolution of the 
width (Full Width Half Maximum) of the same ^^Co lines. Because the Comp- 
ton cross section decreases slowly with energy and expansion effects dominate, 
the evolution of the line displacement and width is relatively insensitive to the 
specific line. Figs. 6 and 7 show that the rate of decline of the blue shift of the 
0.84 MeV ^^Co line is only slightly faster than that of the the 1.24 MeV line and 
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Figure 7: Time evolution of the Full Width Half Maximum (FWHM) of the 0.84 

and 1.24 MeV Co lines. Note that the FWHM is insensitive to the specific line 
considered. 



that the evolution of the line widths are virtually identical for all models. This 
common behavior might make it easier to derive constraints from low signal to 
noise data from missions like CGRO and INTEGRAL. The line shift is critical 
for a proper analysis of proposed high sensitivity but narrow wavelength band 
instruments such as Laue telescopes. 

During the first few weeks, the maxima of all 7-ray lines are blue shifted by a 
velocity that is typically 5000 to 8000 km sec"-*^ with respect to the rest wave- 
length because the ejecta are optically thick to 7-rays during this epoch and the 
flux arises from the approaching side. Photons escape essentially radially from 
the ^"Ni region that is expanding towards the observer. The blue-shift is compa- 
rable to the intrinsic line width, another measure of the kinematic velocity of the 
radioactive matter. The rate of decrease of the blue-shift depends on the optical 



depth of the "^"Ni region and of the matter above it. Consequently, the displace- 
ment tends to drop most rapidly with time in the sub-Chandrasekhar models 
HeD6 and HeDlO, and most slowly in the envelope models such as DET2env4 
or the pulsating delayed detonations such as PDD6 (cf Fig. 6). The evolution 
of the blue shift during the early epochs, prior to 30 days after the explosion, 
roughly two weeks after optical maximum, might allow discrimination of Chan- 
drasekhar and sub-Chandrasekhar mass models of the subluminous variety of 
SNc la as can be seen by comparing the predicted displacement evolution for 
models PDD5 and HeD6 in Fig. 6. With the exception of early time observa- 
tions, the line shift is not a powerful discriminant between helium detonation and 
delayed detonation models. Careful observations of the blue-shift over the first 
150 days may, however, provide such a test for models with masses larger than 
Mq^. If observed, a slow evolution of the blue shift could prove the existence of 
the shell-like structures expected for mergers. 

Figure 7 shows that 7-ray line widths can provide detailed discrimination of 
models for instruments with spectral resolution of order of X/dX ~ 200 — 300. 
Since the line width depends on the velocity distribution of the radioactive ma- 
terial but is insensitive to the wavelength, the results for the 0.84 and 1.24 MeV 
lines of ^^Co given in Fig. 7 are quantitatively valid for other lines as well. As 
shown in Fig. 7, the line width for most models increases on time scales of 10 to 
100 days because radioactive material from an increasing velocity range becomes 
visible as the 7-ray optical depth declines. Since more Ni generally implies 
less shielding by overlying layers and, thus, higher escape probability at early 
times within each series of models, the rate of increase of line width is smallest 
for models with lower ^^Ni masses (e.g. DD201c vs. DD202c). For all models, 
the rate of increase in line width is also smallest for models with "well hidden" 
^^Ni . This is the case if the model naturally produces a definite shell-like struc- 
ture such as PDD6, or if the total mass is larger than M(jj^ as for the merger 
models like DET2env4. 

The behavior of the sub-Chandrasekhar mass models with detonating outer 
helium shells is also conspicuously different in terms of the line width behavior. 
The lines are very broad beginning soon after the explosion, because of the high 
expansion velocity of the outer ^"Ni produced by the helium detonation. The 
line width decreases rapidly at about 2 weeks after the explosion. The rapid 
decline in line width in the helium-detonation models can be understood in terms 
of the line profiles illustrated in Figs. 4 and 5 as being due to the increased 
emission from the narrow inner ^^Ni component (see above). With increasing 
time, this narrow component from the centrally-produced Ni dominates. The 
secondary minimum occurs approximately when the monochromatic flux due 
to the central ^^Ni at a low projected velocity dominates the flux of the high 
velocity component shifted by one Doppler width as the contribution of the high 
velocity component becomes insignificant due to the larger velocity spread and 



smaller mass. For normally bright helium-detonation models such as HeDlO, the 
central Ni mass exceeds that of the outer region by more than a factor of ~ 7 
whereas for subluminous models such as HeD6, this factor is less than 3 (Table 
1). Consequently, for subluminous helium detonations, the early spike in the line 
width is especially pronounced and should be a clear diagnostic for such models. 

After about 50 to 150 days, depending on the model, the line widths approach 
asymptotic values as the envelopes become transparent to 7-rays. In general, the 
asymptotic value of the FWHM line width corresponds neither to the expansion 
velocity of the matter within which Ni dominates the abundances nor to the 
mean velocity of the Ni . Rather, it measures the projected expansion velocity 
at which half the amount of radioactive material is located relative to the material 
seen with zero expansion velocity. For models with centrally concentrated ^^Ni , 
the FWHM depends mainly on the amount of ^^Ni and the size of the central 
hole of the ^^Ni distribution. The former can be independently determined by 
the line flux at late times (see next section). Therefore, the line width provides 
a unique opportunity to probe the very central region and to get an insight into 
the initial conditions at the time of the explosion (see above) , namely, the density 
where the burning starts and the central density of the white dwarf. 
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Figure 10: Same as Fig. 8, for ^^Ni 
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5. Gamma-ray Line Fluxes and Light Curves 
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Figure 8: Time evolution of the integrated line fluxes of the most prominent 

gamma-ray lines of ^^Ni and ^^Co as predicted by theoretical models (see Table 
1) assuming a distance of lOMpc 



The 7-ray line fluxes and resulting light curves provide another diagnostic to 
choose among models for the explosion. The information available will depend 
on whether one knows, from other means, the distance to the source. Fig. 8 
presents the line flux light curves for the two strongest lines each of Ni and 
^^Co for a range of relevant models. Fig. 9 gives the same information for ^^Co , 
including some of the weaker lines and Fig. 10 does the same for ^^Ni . 

We first consider the information contained in distance-independent quanti- 
ties, such as flux ratios. Fig. 8 shows that the helium-detonation models with 
radioactive material close to the surface are predicted to have an especially strong 
line flux at early times. For helium-detonation models, the maximum absolute 
flux due to the Ni lines exceeds that of the strongest Co lines by about a 
factor of 5. For all other models, the Ni flux remains comparable to or smaller 
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Figure 9: Same as Fig. 8, for Co, but including wealcer lines. 



than the flux in ^"Co lines. This provides a crucial, distance- independent test 
for the nature of SNe la. 

Line flux ratios for lines of the same radioactive isotope provide a direct mea- 
sure of the Compton scattering beyond the region containing the radioactive 
species and thus another way to distinguish different scenarios and to separate 
explosions of sub-Chandrasekhar and Chandrasekhar mass progenitors. This 
is illustrated in Figs. 9 and 10. Depending on energy and model, the fluxes of 
^^Ni lines peak between 10 and 20 days after explosion and those from ^^Co peak 
between 50 and 90 days. 

Figures 3 and 9 show that near and after optical maximum, the ^^Ni lines at 
energies below 500 keV have strengths comparable to the strongest 7-ray lines at 
.75 and .81 MeV. The positron annihilation line at .511 MeV is only moderately 
weaker than other lines in the MeV range. These results raise the possibility. 



depending on the sensitivity curve of the instrument, that low energy ^"Ni and 
the 0.511 MeV hues may be detectable more easily than lines in the MeV range. 
At early times, all positrons are trapped (Colgate et al. 1997) and so the results 
do not depend on model-dependent assumptions that would alter the results of 
positron trapping at later times. 
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Figure 11: Luminosity above lOkeV relative to the energy production by radioac- 
tive decays. 



Measurement of line fluxes at both early and late times provides another 
distance- independent means to test models. Fig. 11 presents the ratio R be- 
tween the total 7-ray luminosity and the total instantaneous energy release by 
radioacive decay. The ratio R measures the transparency of the expanding mat- 
ter above the region of radioactive elements. It is very model dependent at early 
times, but eventually approaches a value of unity 4 to 6 months after explo- 
sion. Models with high masses and shell-like structures show small R at early 
times whereas the helium-detonation models with masses less than show 
the highest line ratios. 

Besides a comparison of the line strengths of ^^Co and ^^Ni , the change of 
the line ratios provides important information on the structure of the envelope. 
This constraint comes by comparing lines at different energies as illustrated in 
Fig. 12 which presents results for the 1.24 and .84 MeV lines of Co for a 
range of models. By comparing critical line ratios with their asymptotic values, 
the Compton depth can be constrained. Fig. 12 shows that the ^^Co line 
ratio is sensitive to both time and the model. The line ratios depend mainly 
on the amount of material above the radioactive layer; however, self-shielding 
within the radioactive region is also important. At early epochs, even the helium- 
detonation models with their exposed outer layers of Ni show line ratios that 
substantially exceed the asymptotic value (Fig. 12). Observations during the 
first 2 to 3 months are most suitable in this regard. To distinguish helium- 
detonation models from the merger scenario and the explosion of massive white 
dwarfs, the line ratios must be known to better than ^ 20 percent. Note that 
late time observations are not actually required because the asymptotic value is 
fixed by the branching ratio (see Table 1). 
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Other constraints on the models are possible in the case where the distance is 
known from, e.g., 6 Cepheids (Saha et al. 1996, Freedman et al. 1996) or from 
optical light curves and spectra of (e.g. Hamuy et al. 1996, HK96, Riess, Press & 
Kirshner 1996). In this case, models can be distinguished by the time evolution 
of absolute line fluxes as illustrated in Figs. 8, 9, and 10. For the helium- 
detonation models, the flux in the strongest Ni lines should be in excess of 
10^ photons sec^-*^ cm^^ if observed at a distance of 10 Mpc as the flux peaks 
between 10 and 20 days after the explosion. Late time measurements would 
allow the determination of the absolute amount of radioactive material for given 
distance. 

5.1. Determination of the Time of the Explosion 

The early, rapid change of the spectra may allow an accurate measure of the 
time of the explosion (Clayton et al., 1969). Photons at similar energies encounter 
the same degradation due to Compton-scattering and, thus, the time-dependent 
flux ratio of the lines of similar energy does not depend on the model as shown 
in Fig. 13. Two spectral regions are especially suitable for determining the time 
of explosion: that at about 0.8 containing the 0.81 MeV line of ^^Ni and the 
^^Co line at 0.84McV, and that near 0.5 MeV containing the positron anhilation 
line at 0. 511 MeV and the ^^Ni line at 0.48 MeV (Figs. 2 and 3). Fig. 13 shows 

r ft 

the ratio of the 0.81 MeV line of '-'"Ni in comparison to the 0.84 MeV line of 
^^Co. The ratio of the 0.48 MeV line to that at 0.511 Mev follows an essentially 
identical trajectory, but with an amplitude decreased by a factor of 3.5. The 



ratio of the ''"Ni to "^"Co lines at 0.8 MeV has the advantage of involving strong 
lines with a large flux ratio, but its use is restricted to times between 11 and 
25 days without being obscured by contributions of Compton scattered photons 
of the Co line. The ratio of the Ni line to the annihilation line at 0.5 Mev 
is based on lines with more moderate fliixes but, in principle, it can be used 
between day 20 and 50 for nearby SNe la. 
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Figure 13: The ratio between the ^^A^i [0.81] /^*^Co[0. 84] lines. This ratio depends 
sensitively on the time after the explosion, but it hardly depends on the model. 

The line ratio between the 0.48 MeV Ni and the 511keV positron annhilation 
lines is 3.5. 



5.2. Effect of Integration Times 

Integration times for 7-ray lines may be on the order of days to weeks. There- 
fore it is useful to characterize the influence of the integration time on the result- 
ing spectrum which depends on both the mean time and the underlying model. 
The model dependence enters because the time evolution of the absolute fluxes 
determines the weighting of the line fluxes over the observed time interval. 

Fig. 14 illustrates the change of the spectrum with integration time for the 
specific example of the delayed detonation model DD201c at about 14 and 23 
days after the explosion. For At smaller than 10 days, the resulting error in 
the line flux remains smaller than 10 percent for the Chandrasekhar-mass and 
merger models. The moderate size of the error is due to the rather small changes 
in the absolute fluxes compared to the time-averaged spectra (Figs. 1- 5). For 



the helium-detonation models, the situation is different. For this class of models, 
the large early fluxes and rapid change with time (e.g. by a factor of 2 between 
day 14 and 23, Fig. 4) demand integration times less than ^ 5 days to achieve the 
same error level of less than 10 percent. This should not pose a serious problem 
because the high 7-ray fluxes expected from the helium-detonation models would 
allow shorter integration times. Note that the error in the time of the explosion 
can be further reduced if information about the time evolution of the absolute 
flux can be determined directly by 7-ray observations or, alternatively, by using 
restrictions on the model derived by the many constraints outlined in previous 
sections. 
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Figure 14: The predicted spectrum as a function of integration times around day 
14 and 23 for the delayed detonation model DD201c. 



6. Supernova Detection Rate 



The largest obstacle to 7-ray astronomy of supernovae is achieving sufficient 
sensitivity. This directly affects the observability of SNe la with the current 
generation of 7-observatories, i.e. CGRO, the upcoming INTEGRAL mission, 
and proposed third generation of instruments including Compton and Laue tele- 
scopes. The sensitivity of the instruments generally increases for intrinsically 
narrower lines. Since the supernova features tend to be narrower at early epochs 
(Fig. 7), early time observations are highly favored for the detection of 7-ray 
lines. 



The detection limits, F^j^^, depend on the mission and instrument. For the in- 
struments on CGRO, i.e. COMPTEL and OSSE, F^^^ = lO"^-^ photons sec~l 
cm at about 1 MeV (Johnson et al. 1993). This hmit wiU unfortunately de- 
grade to about lO"^'-'- photons sec"-*- cm~^ after the orbit boost that is planned 
as this paper is written. According to its design parameters, the goal of INTE- 
GRAL was to reach (-P^n^^) = 10~^"^ photons sec~-'^cm~ (Von Ballmoos et 
al. 1995). Recent tests have, however, suggested that this limit must be revised 
upward. Assuming an integration time of 10^ seconds and a 3 u detection, the 
line sensitivities are 10~^-' and 10~^-^ photons sec cm for a line width 
5E /E of 1 and 10 % , respectively. For the classical 7-ray lines of ^^Co at 0.84 
and 1.24 MeV, the sensitivity will only be comparable to COMPTEL (Mattcson, 
1997). Sensitivities larger by a factor of 2 to 3 can be achieved at energies of 
about 100 to about 700 keV. The continuum sensitivity ranges from about 10~^ 
photons sec"-*^ cm~^ MeV"-*^ at 100 keV to lO"^""^ photons sec"-*^ cm~^ 
MeV~l at about IMeV. 

The third generation of proposed 7-ray instruments are likely to use either the 
Compton effect or Laue refraction. The first type of instrument reproduces the 
well-proven technology of COMPTEL. The second type uses the new technology 
of focussing lenses. The advantage of the first type is a wide energy range that 
makes it the ideal multi-purpose instrument. The disadvantage is its size which 
boosts problems both with the background reduction and with the price tag. 
Laue telescopes can, in principle, overcome both of these disadvantages because 
of their compactness. Such instruments are, however, only able to cover a very 
narrow energy range comparable to the intrinsic line widths predicted for super- 
novae. Thus they can detect only one line feature at a time. Currently, the Laue 
instruments must be regarded as an option to study individual selected lines. For 
both types of instruments, the realistic design goals are a broad line sensitivity 
of the order of 10~^ photons sec"-*^ cm~^ (Kurfess 1997), substantially more 
sensitive than CGRO or INTEGRAL. 

The observability of supernovae depends on both the rate of of occurence of 
events and their distance. In recent years, the number of supernovae discovered 
nearby has been rather constant, suggesting that essentially all nearby super- 
novae that could be targets for 7-ray spectroscopy are being discovered. Table 
2 gives the number of SNe la discovered in 1991 and 1992 based on the Asi- 
ago catalog. Only a small fraction of these events have been analyzed in detail; 
however, most SNe la are of similar brightness, comparable to normally bright 
SNe la. From recent studies, the mean absolute magnitude at maximum in V is 
^ -19.3"^ with a spread of about 0.4"^ (e.g. Hamuy et al. 1996, HK96). Sta- 
tistically, the absolute brightness can be translated into an estimate of distances 
to individual observed SNe la as presented in Table 2. 



From the general properties of observed supernovae and the light curves and 



TABLE 2 



Observed SNe la from the Asiago catalog and the expected 7-ray flux in the 0.84 MeV hne of 
^^Co based on the delayed detonation model DD201c. To get estimates for other models the 
results can be scaled according to Fig. 8. Maximum flux in the model is at 70*^. 
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spectra of individual events, the delayed detonation and pulsating delayed deto- 
nation models seem to be the most promising ones to account for the majority 
of SNe la events (e.g. Hoflich et al. 1997 and references therein). Therefore, 
we have chosen the delayed detonation model DD201c to provide a guide line 
for the observability of SNe la in 7-ray lines. For all models of normally bright 
events, the maximum 7-ray flux varies only within 50 percent (see Fig. 8). 
Since the flux scales inversely with the distance that effect will dominate, and 
our particular choice of model does not substantially effect this discussion. For 
SNe la at the luminous end of the observed luminosity distribution, distances up 
to which a detection may be possible should be scaled up by, at most, 50 percent. 

For CGRO after the orbit boost with a detectability limit of F^j^^ ^ 10~^'^ 
photons sec~^ cm— 2, Table 2 shows that normally bright SNe la will be de- 
tectable up to a distance of 7 Mpc with an upper limit of about 10 Mpc. Statis- 
tically, there is an SNe la at a distance between 4 and 9 Mpc about once every 
8 years (e.g. 4.6 Mpc for SN1972e, 4.5 Mpc for SN1986g, 8.4 Mpc for SN1989b) 
somewhat higher than the simple estimate of Table 2. Such an event would 
clearly allow for a positive measurement of the 7-ray fluxes or would rule out 
all the present models. Due to the long expected operational time of CGRO (> 



10 years), there is a good chance that a clear detection of ''"Co lines could be 
obtained which would allow different models within each explosion scenario to be 
distinguished. The low energy resolution of CGRO would not permit the use of 
information on the velocity distribution. If, however, the relative fluxes of, e.g., 
the 0.84 MeV and 1.24 MeV lines could be measured to within an accuracy of 
10 percent and observations are available up to 50 days after the explosion, the 
line ratios would provide a distinction between the explosion of a massive white 
dwarf and that of a sub-Chandrasekhar, helium-detonation. 

A much cleaner distinction between Chandrasekhar and sub-Chandrasekhar 
models could be achieved if observations were available at early epochs, about 
15 to 20 days after the explosion corresponding to about maximum visual light. 
Several SNe la have been discovered up to about 2 weeks before maximum (see 
e.g. Hamuy et al. 1995 for a list), so 7-ray observations at these times may not be 
unrealistic. Figure 10 shows that the peak ^^Ni flux of the sub-Chandrasekhar 
models exceed that of the Chandrasekhar mass models by a factor of 5 for nor- 
mally bright models about 10 days after the explosion and by factors of several 
at later epochs. At 10 days after the explosion, the normally bright helium det- 
onation models predict a fliix of ~ 2x10 photons sec cm at 10 Mpc, 
and at maximum light, 15 days after the explosion, a flux of ^ 10~^ photons 
sec cm ^, factors of about 5 and 2, respectively, brighter than predicted at 
70 days. Given the high maximum ^^Ni line fluxes of helium-detonation models, 
both normally bright and subluminous Type la should be detectable by CGRO 
up to distances of 15 and 10 Mpc, respectively, if the observations start about 
1 week before maximum light in the optical wavelength range or these scenarios 
can be ruled out. SNe la within 15 Mpc occur about once every 5-10 years. In 
particular, the helium-detonation scenario predicts the detectability by CGRO of 
all normally bright SNe la before optical maximum in the Virgo cluster at about 
18-20 Mpc where an SNe la should have an apparent magnitude of V=12-12.5. 
Based on Table 2, the opportunity for such an observation might be available 
every 3 to 4 years as supernova searches of all kinds become more successful at 
discovering SNe la well before maximum light. In particular, SN 1991T would 
have provided a rigorous test of the helium-detonation models had the CGRO 
observations been made at a sufficiently early epoch. The key point is thus if 
an ppropriately early CGRO observation significantly fails to detect a flux of 
^ 5x10 photons sec cm from an SNe la of V=12.5, then helium det- 
onation models can be ruled out. A detection or even a flrm upper limit for 
any SNe la in the Virgo Cluster will provide a much needed distinction between 
Chandrasekhar and sub-Chandrasekhar models. 

With INTEGRAL, the next generation 7-ray observatory, all the types of ob- 
servations discussed here are feasible for the classical strong ^^Co lines which 
have been discussed in the context of CGRO. In addition, the high velocity reso- 
lution (fti 600 km sec~^ , Matteson 1997) will allow detailed studies of line proflles 



(see above). An especially exciting possibility arises from the high sensitivity at 
low energies and the low energy threshold of INTEGRAL which will allow de- 
tections of both the 0.511 MeV line of positronium and the ^^Ni line at 0.48 
MeV. Although the corresponding lines are weaker by a factor of 4 compared 
to the strongest lines of ^^Ni and ^^Co , the increased low energy sensitivity of 
INTEGRAL should allow detection of the positron decay line up to about 10 to 
13 Mpc for normally bright SNe la over a period of several months. Observations 
during the first month after explosion would allow detailed studies of Ni lines 
and the 0.511 MeV line and, consequently, would provide for the first time a re- 
alistic way to determine the time of the explosion (Fig. 13). With INTEGRAL, 
the very strong ^^Ni lines at about 0.16 and 0.27 and the Compton continuum 
can be detected (Fig. 1) out to distances including the near side of the Virgo 
cluster. Therefore, a target of opportunity could be expected every 3 to 4 years. 

For the third generation instruments, all SNe la within a distance of 70 to 100 
Mpc would be within reach. This translates into about 6 targets of opportunity 
per year. Detailed studies of a variety of supernovae including their statistical 
properties would be possible. All the effects mentioned in the previous sections 
could be studied in great detail, opening a new window to do 7-ray astronomy 
of supernovae on a regular basis. 

The sensitivity of the proposed third-generation experiments brings an entirely 
new aspect to 7-ray astronomy of supernovae. Up to this point, the current dis- 
cussion has been based on the premise that all supernovae to be studied will be 
discovered by standard ground-based optical search techniques. The record of 
distant supernovae discovered in this way must be very incomplete. As can be 
seen from Table 2, the number of discoveries does not scale with the volume, 
even for distances not effected by clustering. This strongly suggests the pres- 
ence of selection effects which must depend on the Type of the host galaxy, the 
extinction, and the distance from the center of the galaxy. These questions are 
highly relevant to clarify the progenitor population and its influence on the rate 
of discovery and nature of the supernovae. These issues have important conse- 
quences for the use of SNe la to obtain extragalactic distances at high red shift 
(Hoflich et al. 1997). 

7. Positron Trapping 

One of the more interesting and still unsolved problems of SNe la is the physics 
that determines the slope of the late-time light curve. As emphasized by Colgate 
(Colgate, Petschek & Kriese 1980, Colgate 1983, Colgate et al. 1997), the expo- 
nential decay continues for many e-folds, in excess of 600 days (Kirshner and Oke 
1972, Kirshner et al 1993). There is no observational evidence for the decline 
of the optical light curve due to an "infrared catastrophe," (Weaver, Axelrod 
& Woosley 1980, Fransson et al 1996) that is predicted to occur as standard 



thermonuclear models expand and cool. With assumed complete trapping of 
positrons due to tangled magnetic fields and the resulting small Larmor radii, 
the line emission shifts to the IR fine structure lines of iron. Colgate has argued 
that this problem requires a complete shift of paradigm and in favor of a model 
based on core-collapse that ejects only a small amount of highly nickel-enriched 
matter (~ OAMq) Colgate proposes that with the small mass and associated 
radial combing of the magnetic field, the positrons would leak out with a depo- 
sition function that is essentially self-similar to that for 7-rays. He argues that 
the resulting deposited power in positrons could be restricted to optical emission 
through fluorescence effects and could then give the observed long time expo- 
nential decay. This model is very unlikely to apply to SNe la since it would give 
spectra anomolously rich in ^^Ni and poor in intermediate mass elements, con- 
tinuua that are too blue (for the same reasons as the helium-detonation models) , 
excessively high velocities, light curves that peak and decline too rapidly, and 
implied distances that are incommensurate with those determined by 6 Cephcids. 
Nevertheless, Colgate is right that the problem of the late-time decay must be 
solved in the context of the thermonuclear models before they can be regarded 
as complete. Gamma-ray observations should help to clarify this important issue 
of physics. 

To constrain positron trapping, late-time observations are necessary. Even 
with no magnetic field, positrons are not predicted to begin leaking, thereby 
reducing the positron energy deposition function, until about 300 days. The flux 
is correspondingly small at this time. INTEGRAL would require a supernova 
within only a few megaparsecs to make the corresponding detection. The next 
generation of 7-ray telescopes could, however, make a substantial contribution to 
this issue. A flux limit of 10 photons sec cm would allow a measurement 
of the difference in the positron flux between fully trapped and freely leaking 
models for supernovae out to the Virgo Cluster. 

The issue of the positron leakage could be approached in several ways, some 
of which might be distance dependent and some of which would require an inde- 
pendent measure of the distance. One method would be to look carefully at the 
0.511 Mev line flux, in an absolute sense, or in comparison to the ^^Co lines at 
the same epoch. Another technique, which would not require narrow line resolu- 
tion would be to look at the absolute or relative photon flux of the positronium 
continuum. Because the photon flux is constant it may even be detectable as a 
constant contribution to the the X-ray fliix by X-ray instruments due to their 
superior sensitivity. 

8. Discussion and Conclusions 

We have presented 7-ray spectra and light curves for a variety of models 
of SNe la including delayed detonation, pulsating delayed detonations, merger 



scenarios and helium detonations. Gamma rays trace the radioactive isotopes 
whereas optical photons trace elements and, therefore, 7 observations must be 
seen as an important complement to those at other wavelengths. 

A detailed discussion of the properties of the emitted spectra were given. These 
properties provide sensitive tools to study details of the explosion models and 
information which is relevant for the data analysis of the 7-ray observations 
themselves. Since the different model scenarios show very different signatures 
in their time evolution, spectral features and line ratios, the observation of 7- 
rays have the capacity to remove any ambiguity about the basic nature of Type 
la supernovae, whether they involve thermonuclear explosion or core collapse 
and in the expected former case, whether they involve Chandrasekhar or sub- 
Chandrasekhar mass white dwarfs. For nearby supernovae, 7-ray spectroscopy 
has the potential to separate effects due to the initial metallicity of the exploding 
white dwarf, mixing, and nuclear burning. 

A comparison of 7-ray fluxes and bolometric (UV/optical/IR) light curves can 
provide a measure of the escape probability of positrons and, in principle, of 
the magnetic field strength and distribution. Gamma-rays provide an unbiased, 
direct tool to test the calibration of SNe la and to overcome systematic effects due 
to reddening and errors in the calibration of secondary distance indicators. Such 
tests are critical for the future use of SNe la as distance indicators to determine 
the deceleration parameter and other cosmologically relevant parameters. 

SNe la are not a strictly homogeneous group of objects as previously believed. 
Therefore, any study of them must include a large number of events. For 7-ray 
astronomy, this implies that a broad-line sensitivity of about 10~^ photons s"-*^ 
cm is needed to bring within reach all SNe la up to a distance of 70 to 100 Mpc. 
Based on current discovery rates in the optical, this translates into about 6-8 SN 
yr"-*^ each of which would be detectable for a period of about 6 months. 

The increased sensitivity of future generations of instruments will give for 
the first time the absolute and the relative rate of SNe la as a function of the 
galaxy type unbiased by extinction and the type of the host galaxy which, as 
can be seen from Table 2, hampers optical searches. A monitoring of distant 
clusters on a regular basis is feasible, because an entire cluster can be observed 
simultaneously given the wide field of view of gamma ray detectors. This should 
give valuable information on the progenitor evolution, still one of the major 
mysteries associated with SNe la. 
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